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Abstract 

Using the instant form dynamics of Poincare invariant quantum mechanics and 
the modified relativistic impulse approximation proposed previously we calculate 
asymptotics of electromagnetic form factors for the deuteron considered as two- 
nucleon system. We show that today experiment on the elastic ed-scattering has 
reached asymptotic regime. The possible range of momentum transfer when the 
quark degrees of freedom could be seen in future JLab experiments is estimated. 
The explicit relation between the behavior of deuteron wave function at r = and 
the form factors asymptotics is obtained. The conditions on wave functions to give 
the asymptotics predicted by QCD and quark counting rules are formulated. 
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Recent experiments with electron accelerators and, especially, the Jefferson Laboratory 
(JLab) experiments using continuous electron beam arouse the interest in the theoretical 
study of the hadron electromagnetic structure (see, e.g., [I1121E1I11E1E1E1IH1E] and 
the review [10]). One of the most impressing JLab recent results is the discovery, in 
the experiments on polarization transfer in electron-proton scattering, of a new ("non- 
Rosenbluth") behavior of proton form factors, that is of the faster decrease of proton 
charge form factor as compared with the magnetic form factor [11]. The other one is 
the measurement of the deuteron tensor polarization component T2o(Q^) (Q^ = — g^, q 
is the transferred momentum) in polarization experiments on elastic electron-deuteron 
scattering at ~ 2 (GeV/c)^ [12j. In [13] we have shown that the existing data for 
T2o(Q^) provide a crucial test for deuteron wave functions. 

The JLab program of investigations on elastic electron-deuteron scattering at ~ 
10 (GeV/c)^ [13] attracts exclusive attention. Those experiments will be an important 
source of information on nucleon-nucleon interaction at short distances, on the role of 
relativistic effects, on exchange meson (two-particle) currents and also on the quark de- 
grees of freedom in the deuteron electromagnetic structure. There exists a hope that 
these JLab experiments will help to determine the limits of the two-nucleon model and 
to clarify the interplay between nucleon-nucleon and quark approaches to the deuteron. 
The momentum transfer range ~ 10 (GeV/c)^ (as we show below) is asymptotical for 
the deuteron considered as a nucleon-nucleon system. That is why these experiments are 
of great interest. In fact, first, in the asymptotical domain there exists the most surely 
established QCD prediction for deuteron form factors [15] (see the quark counting rule 
prediction, too [16j). Second, the asymptotic behavior of the nucleon phenomenological 
model probably can give a possibility of correct taking into account of quark degrees of 
freedom in the framework of nucleon-nucleon dynamics at short distances. 

The present paper is devoted to a theoretical investigation of asymptotic behavior of 
the deuteron form factors at large momentum transfer in the framework of the nucleon 
model of deuteron. We use the relativistic invariant impulse approximation in a variant 
of instant form dynamics of Poincare invariant quantum mechanics (PIQM) developed in 
our papers previously ^TTJ dHl dSl 120] and the results of our paper [21] where a theorem on 
asymptotic estimation of the quantities in question is proved which is valid in relativistic 
as well as in nonrelativistic case. 

Let us consider first the asymptotics of the deuteron form factors in the nonrelativis- 
tic impulse approximation. The standard expressions for electromagnetic deuteron form 
factors can be written in terms of wave functions in momentum representation in the 
following way (see, e.g. [22]): 

G'^''(Q') = E / dkk'^ dk' ui{k)gc{k ,Q\k')uv{k') , 
1,1' 

GViQ") = -^T. k'dkk"dk'u,ik)g^^Qik,Q\k')u,ik'), 
^ 1,1' 

GZ''{Q') = -MdYl I kUkk"dk'ui{k)~gl^lj{k,Q\k')uiik') . (1) 
1,1' •' 

Here Gf^ , i = C,Q, M are charge, quadrupole and magnetic dipole nonrelativistic form 
factors, ui{k) are model deuteron wave functions in momentum representation, /, I' = 0,2 
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are orbital angular momenta, Md is the deuteron mass, g^l 



C,Q,M are nonrela- 



tivistic free two-particle form factors (2x2 matrices). The explicit form of g^^, which is 



rather cumbersome, is given in the Appendix of 

In pT] the authors, on the base of a theorem proven there, obtained for the form 
factors ([T]) the asymptotic expansions in inverse powers of Q"^ at Q"^ ^ oo. Two leading 
terms have the form: 



aV2 d 
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(2) 



t=0 
t'=0 



with 



l,l'=0,2 

i = C,Q,M, Ac = 1, Aq = 2MJ/Q^ Am = -Ma, the function t'{t) describes the 
boundary of the domain of integration in ([T]) (see |13j), the parameter tq is defined by the 
deuteron matter radius in different models of NN- interaction |21j . 

In modern calculations the deuteron wave functions usually are of the following ana- 
lytic form in the momentum representation (see, e.g. 



M^^) = /;2 , 2A' 

V TT ^ (fc^ + mj) 



U2{k) 



or in the coordinate representation: 

Uo{r) = ^ CjCxp {—rrij r) , U2{r) = ^ D^exp {—rrij r) 



1 + 



+ 



nij r {rrij ry 



(3) 



(4) 



here rrij = a + mo (j — 1) , a = y M \ed\, M is average nucleon mass, Ed - the binding 
energy of the deuteron. The coefficients Cj, Dj, the maximal value of the index j and mo 
are determined by the best fit of the corresponding solution of the Schrodinger equation. 
The standard behavior of the functions at short distances: 



Uo{r) ~ r , U2{r) ~ , 
is provided by imposing the following conditions on the coefficients in (II]): 



(5) 



(6) 



Using the explicit form of g^- 



we obtain from 
factors: 



r^NR 



1 2 
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C, Q, M from [13] and taking into account the Eq. 



E^X 



the main asymptotic terms of the nonrelativistic deuteron form 

2 

{GUQ') + GUQ')) , (7) 
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G 



NR 



1 24^/2 



4 



M 



1 2i5Mrf 



n 2 



m ■ 



(9) 



here G^^^l^-^ are charge and magnetic form factors of proton and neutron, correspondingly. 

Let us note that the asymptotic forms ©-Q are vahd for any model of A^A^-interaction 
if (E]) is satisfied. The fact that Gq^ decreases faster than other form factors is due to 
the faster decreasing of the D-wave function at r — >■ as compared with S'-wave function 
([5]). From the mathematical point of view the form of the leading terms in ©-([H]) is the 
consequence of the conditions (E]). Any modification of these conditions (or, equivalently, 
of the conditions (^) changes ([7])-([n]) as well. 

The asymptotic expansions of the deuteron form factors contain the nucleon form 
factors. We use the standard dipole fit for magnetic nucleon form factors [TOj : 



Hp 



/in 



1 + 



(10) 



with =0.71 GeV^. Here fip, fin ^ are proton and neutron magnetic moments, core- 
spondingly. As the neutron charge form factor is small enough we neglect its contribution. 
In this connection, one has to notice that its behavior at large momentum transfer has 
also the dipole form (see, e.g., |12]), so that its contribution would not change our results. 
We assume the " non-Rosenbluth" [H] behavior for the proton charge form factor and so 
we obtain the faster decrease of the deuteron charge form factor as compared with the 
magnetic one. Let us note that the asymptotics of nucleon form factors (|TU|) (~ Q^'^) 
agrees with the predictions of QCD for nucleons [T5] . 

Usually describing the deuteron one is dealing with the asymptotic behavior of the 
function (see, e.g. .[14]): 

FAQ') ^ 



, (11) 

where the structure function A{Q'^), which enters the electron-deuteron scattering cross- 
section, can be written in terms of the deuteron form factors in the following way [lOj : 



Gim 



Iv'Gim+lvGiAQ'' 



V 



4MJ 



(12) 



After substituting the expansions ([7])-([nD in (|T^ one can see that the main contribution 
to (fTTl) or (fT2l) comes from the term containing the deuteron magnetic form factor. Note, 
that the "non-Rosenbluth" behavior of proton form factor reinforces the main role of 
this term. Using (ITOl) we obtain the asymptotics of AIQ"^) in the nonrelativistic impulse 
approximation: 



1 2^^{fip + fin) 
37rr|]M2 



E Gjm^ 



(13) 
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Let us emphasize that our approach permits us to write the prefactor in (fT3ll exphcitly. 
Other authors deahng with the deuteron asymptotics (see, e.g., [26J), usually derive the 
power of momentum transfer only. Our form of the power-law dependence is 

F^m - Q-'' (14) 

in nonrelativistic case. 

However, at large momentum transfer it is important to use relativistic approach. 
To give the relativistic description of the deuteron we use here the formalism of the 
Poincare invariant quantum mechanics (PIQM). Today this approach is one of the effective 
approaches to the structure of composite systems (see, e.g.,[TD]). Actually we use a variant 
of the instant form dynamics of PIQM developed by the authors (see [IHl UHl [20])) -In 
our approach the relativistic electromagnetic form factors of a composite system can be 
written in a form similar to the nonrelativistic case ([T]). The equations obtained in the 
modified Lorentz-invariant relativistic impulse approximation are given in the paper [13] 
and have the following form: 

G^iQ') = J2[ dVrsdV^ipiis)gj^cis,Q\s') Ms'), 
1,1' 

9 r 

^ 1,1' 

G'^iQ') = -M,J2 I d^sd^'Ms)g'^^{s,Q\s')Ms')- (15) 
1,1' 

Here ^i{s), 1,1' = 0,2 are the deuteron wave functions in the sense of PIQM, g^- , i = 
C, Q, M are relativistic free two-particles charge, quadrupole and magnetic form factors 

m- 

The deuteron wave functions in the sense of PIQM are solutions of the eigenvalue 
problem for a mass squared operator for the deuteron M| \ip) = Mj \ip) (see, e.g., [TU|[TH]). 
This equation coincides with the nonrelativistic Schrodinger equation within a second 
order in deuteron binding energy e^/4M (M is an averaged nucleon mass). The value of 
this quantity is small, so the deuteron wave functions in the sense of PIQM differ from 
nonrelativistic wave functions by the normalization only. In the relativistic case the wave 
functions are normalized with relativistic density of states: 

roc rjU 

E^y^ ^f(fc) ^^7=^ = 1, y^i{k)= \rsku,{k), s = 4{e + M'). (16) 

The nonrelativistic formulae ([T]) given in [22] , can be obtained from relativistic ones ( ITSj) 
in the nonrelativistic limit. 

Two leading terms of the asymptotic expansion of the deuteron form factors have the 
form [21j (see also the notes to Eq. ((21)): 

GfiQ') - A^^^^ f5//^(t,Q^t'(t)) + ^|;//^(t,Q^^'(t)) + 
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8M ^2 ^r., ^2 // 



t=0 
t'=0 



(17) 



where 



LV=o,2 4vss 



V 2 V 2 

Using ([3]), ( fT6l) . ( JTTl) and the exphcit forms of g^l , i = C,Q, M from [13], we obtain 
the main terms of relativistic asymptotic expansions f|T5l) . f|T7|) . which enables us to write 
the relativistic asymptotics of the function AIQ"^) in ffT^ in the form: 

A'iQ') - ^.A^^'iQ')- (18) 

So, the relativistic asymptotics of the form factor (fTTi) is: 

F^iQ') ~ Q"^ . (19) 

The relativistic effects slow down the decreasing of the form factors as compared with the 
nonrelativistic case (fT3l) . ( |T4|) . 

The comparison of the first and the second leading terms in the obtained asymptotic 
expansions for different models [231 [2H of A^A^-interaction shows that for all of 

them the asymptotic form (|T^ . (HM is valid at — 6 (GeV/c)^. So, the existing JLab 
experiments at large momentum transfer [28] have already reached the asymptotic region 
predicted by relativistic two-nucleon deuteron model. 

In this connection it is of interest to compare directly the obtained asymptotic predic- 
tions with experimental data. Fitting the existing experimental points for seven highest 
attained values of momentum transfer [28.| in region 3.040-5.955 (GeV/c)^ by a power-law 
function we obtain the following estimation for ( fTTi) with = 3.93 ■ 10~^: 

So, comparing (fT^ . (fT^ and (^01) one can see that up to fitting accuracy, the experimental 
data are described by the relativistic formula (|T^ . Let us note that this result does not 
depend on the actual model of A^A^-interaction and in fact is due to general conditions 
([5]) and ([6]) only. So, in the recent JLab experiments the range of momentum transfers, 
which can be characterized as asymptotic one for our relativistic two-nucleon model of 
deuteron, is reached. Let us emphasize that this concerns our relativistic approach only. 
For example, in the relativistic approach of [2S] the asymptotic decrease is faster than the 
experimental one and even faster than that of the quark counting prediction. 

The numerical comparison of our asymptotic prediction ( ITSl) (including prefactor) for 
different model interactions [23l EH |25], [27] with experimental values of ^(Q^) demon- 
strates that our theoretical curves lay lower than the experimental. The same is true 
about our full (without asymptotic expansion) relativistic calculation using Eqs. ([H]), 
( fT2l) . (ITSl) . ( fT6l) . It is worth to notice that the relativistic effects shift the curve closer to 
the experimental one as compared with nonrelativistic case. 
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Figure 1: Experimental data for the deuteron form factor (fTTl) . Dark dots are experimental 
dots obtained in JLab [14J, light dots are projected in future experiment dots [H], solid 
line is our fit of existing experimental data at high momentum transfer (l20ll 



The asymptotics ffT^ . ffT^ are obtained in the framework of the nonrelativistic impulse 
approximation. The asymptotics ffT^ . ffTI?|) correspond to the Lorentz- invariant modified 
impulse approximation [18]. In this connection it is interesting to estimate the asymp- 
totic behavior of the contributions of terms out of impulse approximation, in particular, 
of meson exchange currents (MEC). The agreement of experimental asymptotics of the 
deuteron form factors (120|) with that obtained in the relativistic impulse approximation 
(1191) means that the MEC contribution at large momentum transfer either has the same 
power law dependence on as in (|T^ or decreases faster. It is possible that this fact is 
due to a fast decreasing of transfer mesons form factors (for example of p7r7-form factor) 
at — > oo. This problem will be considered in detail elsewhere. 

Let us compare our results with those of the quark approach and of QCD. At — > oo 
there exists the well established prediction in the framework of these approaches [T^ [TB]: 

FdiQ') ~ Q-'' . (21) 

As one can see, this prediction does not agree with the current experiment (1201) . In other 
words the quark degrees of freedom are not seen in the range of today momentum transfer. 
The hope to see quark degrees of freedom is connected with the possibility of obtaining 
in JLab future experiments the deviation of power law dependence of A(Q^) and Frf(Q^) 
from that given by our relativistic model ffTSl) . ffTI?|) and by today experiment (I2U1) . 

In figure [U the projected results of future JLab experiments are presented [H]. The 
saturation is interpreted [2] as due to quark degrees of freedom. One can see that this 
fact may take place in the range of ~ 10 (GeV/c)^, where the projected data leaves 
our fit (I2U1). 

Let us discuss one more problem connected with the deuteron form factor asymp- 
totics, namely, the possibility to incorporate the QCD predictions in the nucleon-nucleon 
dynamics. One of the possibilities of obtaining the asymptotics ( |2T1) in nucleon physics is 
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to put some conditions on nucleon form factors in (I71)-(l9]). However, this way can not be 
considered as self consistent. More correct way is to use the QCD prediction (~ Q^"^) co- 
inciding at ^ oo with the result of dipole fit ffTUj) used while obtaining the asymptotic 
estimations (fT4l). (fT9|) . 

So, we formulate the problem in the following way: what kind of behavior at small 
distances would have the deuteron wave function, or, in other words, how the nucleon- 
nucleon potential at short distances has to be modified - in order to obtain the asymptotics 
of electromagnetic deuteron form factors predicted by QCD? The answer can be obtained 
easily and the analyzis shows that the quark model asymptotics could be derived in the 
nucleon dynamics formalism if in addition to the conditions ([6]) for the s-wave function 
([3]) the following condition is imposed: 

ECX = 0- (22) 

j 

This condition means that in the vicinity of zero the wave function has the following 
form: 

Uo{r) ^r + ar\ UoiO) = . (23) 

So, we have solved some kind of inverse problem: we found a condition for the deuteron 
wave function to give the asymptotic behavior of the deuteron electromagnetic form factors 
(in the framework of the nucleon model of deuteron) given by quark approach. The 
presence of the quark degrees of freedom changes the deuteron wave function ([3]) following 
the conditions (1^. 

To conclude, in the present paper the following results are obtained. 

The asymptotics of the deuteron form factors at large momentum transfer is derived 
in nonrelativistic as well as in relativistic case in the framework of impulse approximation 
in the PIQM instant form dynamics. 

The relativistic effects slow down the asymptotic decreasing of form factors and result 
in the power-law dependence on which coincides with the experimental data. 

It is shown that the range of momentum transfers of recent JLab experiments may be 
considered as asymptotic for relativistic two-nucleon deuteron model. 

By comparison of the calculated asymptotics with experimental data it is established 
that quark degrees of freedom are not seen in today experiment and could be seen in 
future JLab experiments when the deviation of experiment from the results obtained in 
this paper would be reached. 

The condition for the deuteron wave function at r = to provide the asymptotic 
behavior of deuteron form factors predicted by QCD is obtained. 

This work was supported in part by Russian Foundation for Basic Researches (grant 
No 07-02-00962). 
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